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Abstract 

We analyze the puzzle of the dimuon CP asymmetry in Eg decays in two Higgs doublet models. 
We show that the flavor changing neutral current (FCNC) induced by the Higgs coupling in a 
type ni two Higgs doublet model provides a solution to the dimuon charge asymmetry puzzle by 
enhancing the absorbtive part of the mixing amplitude T\2- We investigate different experimental 
constraints and show that it is possible to enhance in order to explain the dimuon asym- 
metry observed by DO. This enhancement requires large Higgs couplings to the first and second 
generations of quarks which may also explain the recent 3.2 a Wjj excess observed by CDF. 



1 



I. INTRODUCTION 



The DO collaboration has measured the like-sign dimuon charge asymmetry in semilep- 
tonic 6-hadron decays A^i. The following result has been reported 

All = -0.00957 ± 0.00251(stat) ± 0.00146(syst). (1) 

The like-sign dimuon charge asymmetry A^i for semileptonic decay of b hadrons is defined 
as 

A" - <^ - (2) 

where Njj"^ and A^"^ are the number of events containing two 6-hadrons that decay semilep- 
tonically two positive or negative muons with the quark subprocesses: b /x~z/X and 
b jj+uX. 

This result indicates a 3.2 cr deviation from the Standard Model (SM) prediction. A 
confirmation of this deviation would provide unambiguous evidence for new physics (NP) at 
low energy with a new source of CP violating phases. It is a common feature for any physics 
beyond the SM to possess additional sources of CP violation besides the SM phase in quark 
mixing matrix. These new phases can induce sizable contributions to direct and indirect 
CP asymmetries in Bd^s decays and thereby resolve the apparent discrepancies between the 
observed results and the SM expectations. 

The charge asymmetry A^i at the Tevatron can be expressed as A^i = (0.506 ±0.043)0^^ -|- 
(0.494 ±0.043)ag; where the charge asymmetry a^^ {q = d, s) for "wrong- charge" semilep- 
tonic i?°-meson decay induced by the oscillation is defined by 

""si p^^o ^ ^+j,x) + r(50 ^ fi-ux) ■ ^ ^ 

This asymmetry can be written as ^ 

I AF 

where the mass and width differences between Bg mass eigenstates are given by 

AMs, = Mb, - Mb, = 2\M1,\, (5) 
AP, = 2|r?2|cos0„ (6) 

and MI2 and P^g the dispersive and the absorptive parts of the mixing amplitudes, 
respectively. The CP violating phase (pg is defined by (pg = SiTg{—Mf 2/^12)- Iii the SM, 
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the arguments of ^12 aligned at the leading order due to the unitarity of the 

quark mixing matrix, and thus, the phase (pg is small irrespective of the individual phases, 
f3d = aTg{-V:,Vcd/V;,Vtd) and (3, = arg(-V;,Vt,/V;*,K,). The ratio of TlJMf^ is roughly 
proportional to ml/mf with a small phase, and as a result, the SM prediction of the charge 
asymmetry is small: af;(SM) ~ -5 x 10~^, a^^(SM) ~ 2 x 10"^ and A^,(SM) ~ -0.00023 
[4 1, which is in clear contradiction with the DO result in eq.([T]). 

In this respect, it is clear that a large new CP violation in Bs mixing is required to 
enhance a^i if 5^ mixing does not contain a new physics contribution. The phase (pg can 
be large in general in many new physics models because the phase alignment between Mfg 
and can be broken if a new particle, such as a supersymmetric particle, propagates in 
the loop diagram, which contributes to the mixing amplitude js]. However, the absorptive 
part is not necessarily modified by the propagation of new particles due to the on-shell 
condition of the intermediate states. The magnitude of Mfg is determined by the mass 
difference AM^^ , and thus eq. (jl]) tells us that a^i has a maximal value if F is dominated 
by the SM tree level contribution [4]. One can easily show that if one simply extrapolates 
eq.(jl]), then one needs a new phase (pg with siiKpg > 1.6 (which is outside of the domain of 
sin 0s) in order to account for the experimental result of [5,] by the DO collaboration. 
Also with a large smcps ~ 1, the decay width difference AF^ = 2IFJ2I cos0s is suppressed 
and becomes inconsistent with the experimental constraints. Therefore, one concludes that 
the sustainability of the DO results of the like-sign dimuon charge asymmetry in semileptonic 
6-hadrons decay would be a clear hint of possible NP that modify the absorptive part of the 
mixing amplitude F^2- 

In this paper we show that the flavor changing neutral current (FCNC) induced by the 
Higgs coupling in two Higgs doublet models (THDM) can provide a solution to the dimuon 
charge asymmetry puzzle by enhancing F*2- The THDM is classified by the selection of 
the Higgs couplings to fermions. We will consider a general type of coupling (so called 
type HI THDM) to obtain an appropriate FCNC source to modify F^2- We enumerate 
the experimental constraints and investigate if there is room to enhance F^2 to achieve the 
DO result of the dimuon asymmetry. We will obtain the operators generated by charged 
Higgs exchange which can enhance F12. To accomplish this, suitably large Higgs couplings 
to the first and second quark generations needed. Such couplings can, in addition, lead to 
an excess of Higgs decays into dijets. Indeed, the CDF collaboration has recently reported 
a 3.2a excess in the 120-160 GeV range in the invariant mass distribution of the dijets in 
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association with a W boson [6]. The excess may be explained by the Higgs decays to dijets 
via the new Higgs couphngs. 

The presence of two Higgs doublets is required in supersymmetric models and in models 
with the left-right gauge symmetry SU{2)l x S'f/(2)/j x U{1)b-l 01- We briefly explore 
if the left-right symmetric model can explain the large dimuon asymmetry consistent with 
other experimental constraints. 



II. DIMUON ASYMMETRY IN TWO HIGGS DOUBLET MODELS 

The two Higgs Doublet Model is an extension of the SM that naturally introduces a new 
source of CP violation and FCNC jsj. In this class of models, the most general renormalizable 
and gauge invariant Yukawa interactions are given by 

-Cy = Y^^q.u'jH^ + Y:;^q,u'jH, + Y^'q.d'jH, 

+ Y!;'q^d]H^ + Y^n.e'^^Hd + Y^'H.e'jH^, (7) 

where the Higgs fields Hu,d have hypercharges Y = ±1/2, and Hu,d are defined as Hu,d = 
—iT2H*^^. The Yukawa couplings Yf and Yj (/ = u, d, e) are 3x3 matrices with a generic 
flavor structure. In this case, the fermion masses are given as 

Mu = Y:^v., + K^^t;^, (8) 

M, = Y^va + t;., (9) 

Me = Y^^d + Y;:^Vu. (10) 

In the basis where the mass matrix is diagonal, FCNC interactions through the neutral 
Higgs {H and A) exchange are generated from the mismatch between the diagonalization 
of the mass matrices and the Yukawa interactions. For instance, couplings like s^b'^^H^ and 
bLS%H^ are induced by {V^YdVfY^qid'^jHd = {Mf^'^/vd - V^Y^Vf tan ^Y^id^jH a, where 
tan/3 = Vu/Vfii and are the diagonalizing matrices of M^. Hereafter, we will work in 
the basis where is diagonal, and we omit V^, in the expressions. Then, by definition, 

These couplings contribute to Bg-Bg mixing and modify the amplitude M{*2 as follows |9[: 

= m.f'' + (Mi^2)2™, (11) 

where (Mf2)^^° is given by 

y23y32* 

(^12)'™ ^ -jL^^B'^l(^b,sn){bnSL)m. (12) 



m 



H 



We note that even in the minimal supersymmetric standard model (MSSM) this type of 
modification of via the neutral Higgs exchange can be obtained for large tan /3 through 
the finite correction of the Yukawa couphngs due to soft SUSY breaking terms. In this case, 
if = (r^^)^'^, the charge asymmetry a^; is given by 

|/'-pl2\SM| 



sin(</)^^ + 2e,), 



''si 



r.|(Mf)SM| 



(13) 



where r, = |1 + Mf^^/Mf^\ and 26^ = arg(l + M^^^/Mf^). Thus, 26^ can be large if 
l^d^^d^*! ^^sO^d^^P*) l^-rge, satisfying the experimental constraint : ~ 1. As 
mentioned above, a large value of 29s is necessary but not sufficient to account for the 
dimuon CP asymmetry A^i. A significant enhancement of preferable. 

In the 2HDM several Ab = 1 effective operators which may modify F^g) can be gener- 
ated. For example, with non- vanishing Y^^"^ and Y^^^ the operator {hLCR){cRSL) is generated 
through the charged Higgs exchange. Let us enumerate the possible Ah = 1 effective oper- 
ators generated via the Higgs exchange, which are suitable for modifying Ff2. 

1. (6lO(m]?sl) = -|(&l7mSl)(m]?7^0 

3. {bRSL){u'fiui) and {sRhLjiu^Ru'j) 

4. {bRSL){TLTR) and {sRhL){TLTR) 

The operator (6L^iij)('W^SL) is generated through the charged Higgs exchange with the 
coefficient: 



E 



v,tv: 



q' s 



mir. sin^ B 

q,q'=u,c,t tl + 



Y' 



qi 



M*^^ COS 13 



Y' 



(14) 
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As we will see below, this is the preferred operator for modifying V^^- 

The operator (feijti^) (li^s^) is generated through the charged Higgs exchange with the 
coefficient: 



^1 jq 



f— f , m1r+ cos^ P 

q,q'=d,s,b H+ 



-Y' 



q3 



-Y' 
^d 



(15) 



q'2 



This operator can modify through interference with W boson exchange. The effect can 
be large if yj^^ and (y^sin^ — yjcos^)22 are sizable. However, one needs fine-tuning to 
obtain a large contribution to Ffg since the strange mass is = (y^cos^ -|- Y^sm/3)22V. 



Such fine tuning may also generate the operator {bRSL){sLSR), which affects Bd 4>K by 
an amount which causes disagreement with the experimental data. 

The operators {pRSLjiu^RU-'j^) and {suhLjiu'^RU-'j^) can be generated through the neutral 
Higgs exchange. We note that both operators are needed to modify since {u'^fjU''^) cannot 
be self-conjugate. 

Finally, the operators {hj^SL){T jTr) and {sjPl){,'Ti^tr) can be generated from Y^, Yp via 
the neutral Higgs exchange. It is remarkable that when Y^^'^Vu generates the muon mass, 
Br(i?s — 7- rr) is not constrained by experiment and can be modified. Note that in MSSM 
(or THDM with Yg' — )■ 0), the following relation is obtained 



Br(i?s — > rr) / nir 



(16) 



Thus, due to the experimental bound on Br(i?s — > /x/i), the quantity Br(i?s — )• rr) cannot 
be large enough to modify F^g- 

In the case of neutral Higgs exchanges, F12 is modified at one-loop level while A6 = 2 
operator to modify is also generated through the neutral Higgs exchanges at tree level. 
Therefore, to obtain a sizable F^g contribution, a large is also generated, which is 
unwanted. 

Next we consider the following experimental constraints: 

1. B ^ M1M2 and B ^ ££ decays. 

2. 6 57. 

3. lifetime ratio tbJtb^- 

The constraints from two body decays into mesons and leptons are studied in jlO|, and 
the allowed operators for A6 = 1 are 

s6cc, s6fr, dhcu. (17) 

In the MSSM, it is difficult to enhance the A6 = 1 operators selectively because the inter- 
actions are related by known (or constrained) coupling constants. 

The non-holomorphic Yukawa coupling can also contribute to the charged Higgs inter- 
action and lead to an important effect on 6 — )■ 57. In a non-SUSY type II 2HDM, the charged 
Higgs has to be heavier than ~ 350 GeV to agree with the experimental measurement of 
h — )■ 57. However, in our case, this constraint can be relaxed by a small contribution. 
Note that Y^^ {Yd"^) is constrained by hu — )■ 5^7 {pL — ?■ sr'~^) contribution. 
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The 6 — )■ S7 process also constrains the large log contribution due to the renor- 
malization group evolution below the W boson mass ll|. The 6 — )■ 37 constraint 
disfavors the neutral Higgs contributions. This is due to the fact that the operators 
(sL&ij)(/L/ij) and {sRhL){fRfL) (/ = c, 6, r) give leading order large log corrections (5C7,8 ~ 
m f / mi,\n M^r / ml) to the b ^ s'y operators, and hence they are dangerous operators. We 
note that the charged Higgs contribution can generate only next to leading order corrections. 

Finally, the lifetime ratio tbJtb^ provides a stringent constraint on the operator con- 
tributing to r^2- The modification of can induce a change (~ O(10)%) in the lifetime of 
Bsi which can be consistent with the large hadronic uncertainty. However, this uncertainty 
ed in the lifetime ratio tbJtb^- The current world average of the experimental 



is cancel 
result is 



12| 

tb./^b, = 0.99 ±0.03. (18) 



Therefore, if is modified, the lifetime of should also modified, which provides a strong 
constraint. Actually, for system, it is known that the bound for 'B>i{Bfi — )■ rr) < 4.1x10"^. 
Consequently, the modification of FJ2 via Bg — > rr has a deficit due to this lifetime ratio 
constraint. The only allowed operator to modify the B^ lifetime seems to be dhcu, which 
can modify T^^. Fortunately, the operator dhcu alone modifies neither Mj^2 

nor F^2 since cu 

is not self-conjugate. 

The above constraints imply that the only possibility left is the charged Higgs exchange 
operator {qLlfibL){u^B_'y^u-'j^) (g = d,s). The coefficients of these operators are obtained, as 
in eq. ffT^ . and they are proportional to Xif,X*^, where = [(M^^^^ cos f3/v — Y^)^V]ij. In 
order to generate the operators sbcc and dbcu, we need Xcs,Xcb and Xud- The condition 
Xcs ~ Xud can make these effective operators comparable to keep the lifetime ratio tbJtb^- 
The condition can be satisfied when Y(^^ ~ ^u'^- However, if Y^^ is sizable, one needs fine 
tuning to obtain the proper up quark mass. This can be relaxed if tan /3 is large. 

The F^2 contribution is estimated as 



Icc, (19) 



where 7cc = — 4m^/m^(l — 2/3(m^/m^)). The contribution to the dispersive part of the 
mixing can be written as 

^sTHDM XIX^X*,X,, /(mf /M^, m^/M^, M^/Mf,) 



(20) 
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where i,j = c,t and f{x,y,z) and S{x) are the Inami-Lim functions. If we choose FJ^^ 
appropriately, the charged Higgs contribution to can vanish, keeping the contribution 
to r*2, (In fact, if l^*-' oc 6ij and rrii ^ Mh, a GIM-like mechanism works for the dispersive 
part of the meson mixing amplitudes). This is because the top-loop can contribute to MI2, 
but not to r^2- Therefore, one can avoid an excessive contribution to and still modify 
r^2 appropriately in this scenario. 

In Fig{Tl we show r™°^/rf^ as a function of y'^, where we assume yj*-^ = for 
simplicity. We choose the charged Higgs mass to be 160 GeV. 

In order to modify TI2, we need whose magnitude depends on the charged Higgs mass 
rriH- As previously mentioned, we need to fine-tune in order to obtain the proper up quark 
mass in order to satisfy the lifetime ratio. The fine-tuning can be relaxed when the charged 
Higgs is light. However, the constraints from 5" — ?■ tu and B — ?■ Dtu should be taken into 
account [13|. In a general 2HDM, the non-holomorphic coupling ir'^H* can contribute to 
the r mass, and thus the Ir'^Hd coupling (which is important to i? — t- tv) may have freedom 
to relax the B ^ tv constraint. 

The sizable Higgs coupling Y^^ can provide an interesting hadron collider signal. The 
coupling can cause a resonant production of the charged and neutral Higgs bosons (which 
mainly contain Hd)- The Higgs bosons decay into two jets via the Y^^ and Y^'^ couplings, 
which are needed to enlarge T\2- Thus, dijet excesses can be observed around the masses of 
the Higgs bosons. If the Higgs boson is light (< 200 GeV), the dijet events are buried under 
the QCD background. However, even if the Higgs boson is light, there is a chance to observe 
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FIG. 3: r™DM/pSM plotted as a function of F^^s fo^ yni ^ y^22 ^ g.S. 

dijet events produced by Higgs decays associated with W /Z gauge bosons. Recently, the 
CDF collaboration has reported an excess of dijets associated with W boson (which decay 
into iv) [6]. The excess in the dijet mass distribution is in the 120-160 GeV range, and it 
can be explained if the Higgs boson mass is about 150-160 GeV (the peak of the dijet mass 
distribution shifts to lower mass due to cuts). We note that the bottom quark mass should 
be generated by q^b'^H* because 6-quark excess is not observed. 

The charged and neutral Higgs bosons associated with the W boson are produced by 
the t-channel exchange of the left-handed quarks through the coupling. The resonant 
(or off-shell) production of the charged/neutral Higgs boson through can also contribute 
to the neutral/charged Higgs boson associated with the W boson. Similar processes via a 
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Higgs doublet are also analyzed in [ij, ll5|. The CDF estimated production rate of Wjj 
excess is about 4 pb. We use MadGraph/MadEvent (version 5) to estimate the WH cross- 
sections in this model. In FigJJl we plot a{pp — )■ WH) as a function of y'^ for mH+ = 160 
GeV (neutral Higgs masses are also 160 GeV). In the figure we use a K factor ~ 1.35. The 
WH production cross-section is dominated by charged Higgs associated production (roughly 
70%). The charged Higgs boson then decays 100% into jets. We find that the excess can be 
explained for Higgs couplings Y^^^ ~ F^^^ ~ 0.5. We also have ZH (includes both charged 
and neutral Higgs) productions in this model, but the cross-section is at least a factor 3 
smaller compared to the WH production. The ZH production cross-section goes further 
down if the neutral Higgs masses are heavier than the charged Higgs masses. 

In FigJSl we plot T™^^/Tff as a function of Y^^^, where we choose Y^^^ = Y^^^ = 0.5 
and mH+ = 160 GeV. Since in this case we need Y^^^ to be ~ 0.1 to produce the desired 
modification. Due to this small Y^^^, we have a very small amount of single b quark present 
in the Wjj signal. As mentioned before, Y^^^ is constrained by 6/? — t- 5^7 operator. In order 
to turn on the Y^^^ term, one should make Y^^ small and the bottom quark mass should 
come from the non-holomorphic term. The absorptive part of the mixing is insensitive to 
y^^^, and thus Y^^^ can be used to adjust the dispersive part of the mixings. 



III. DIMUON ASYMMETRY IN LEFT-RIGHT SYMMETRIC GAUGE THEORY 



In SUSY models, the non-holomorphic terms [Y^, Y^ couplings) can arise only from the 
finite corrections. A sizable contribution to is easily obtained, but in order to obtain 
sizable effects for F^2 consistence with other experimental constraints such as 6 — 57, O(IOO) 
TeV scale physics should be considered 

In the left-right symmetric gauge theory (LR model), the gauge symmetry is extended to 
SU{2)l X SU{2)ff X U{1)b-l |7|. In this model, a number of Higgs fields and couplings are 
usually needed in order to obtain the required quark masses and CKM mixings. The Higgs 
couplings can generate FCNC, which can affect the dimuon asymmetry, as described in the 
previous section. In addition, the SU{2)ji gauge bosons can also contribute to FJ2 ^^id 



16|. Let us estimate the contributions of these gauge bosons to the mixing amplitudes. 
The RRRR operator (gR7;x6_R) (g^7^&_R) can contribute to F^2 {q = d, s): 

ri2 ^.(Vt,K,\ fg,,\ (M^^\ (l + 0{mllml)\ (21) 



F?f \Vii*Vi \gj \Mw, 
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where and are the left- and right-handed quark mixing matrices, gL,R are the gauge 
couphng constants, and Mw^ r are the masses of the Wl,r gauge bosons (in general, these 
bosons can mix). We have used here the unitarity of the quark mixing matrices: V*gVub + 
VcqVcb + V^*Vtb = 0. Note that the M12 contribution from the RRRR operator is tiny because 
of the unitarity of Vr and mu,c,tl My/^ <^ 1- 

The LLRR operator (gL7At&L)(^7^&_R) can contribute to M^2- 



q 

12 



where the loop functions A2 and F can be found in 17|. Note that the contribution to F 
via the LLRR operator can be negligible because charm or up quark masses are inserted (or 
Wl-Wr mixing is inserted twice). 



As a result, we obtain 



(23) 



The correction to should be less than about 30% when iMf^^/Mf^'^l < 2 and 
M^^/M^^ < 10^'^. Therefore, the contribution via the right-handed gauge boson is not 
a better choice to modify since modification to will be too large irrespective of the 
choice of the right-handed quark mixing matrix. 

The merit of the right-handed gauge boson is to modify T^-^ to adjust the lifetime ratio 
for Bd and Bg. If the right-handed quark mixing matrix is 

^^1 0^ 
1 



V 



R 



(24) 



the operator {d^^ fjhRjicR'y^UR) is generated by Wr exchange, and the lifetime of Bd can be 
tuned depending on gj^/M^^. In the previous section, we found that the element Y^^^ can 
modify Tf^ but it requires fine-tuning to obtain the proper up quark mass. The contribution 
from the right-handed gauge boson can relax the fine-tuning. We must adjust the contribu- 
tion to make the lifetime ratio lie within the current experimental uncertainty in eq.f llSp . If 
we choose the elements (e.g. V^) to be exactly zero, any unwanted contribution to meson 
mixings from the Wr gauge boson can be avoided. 
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IV. CONCLUSION 



In this paper we have investigated the dimuon CP asymmetry in Bg decays in two Higgs 
doublet models. We find that the fiavor changing neutral current (FCNC) induced by 
the Higgs couplings in type III two Higgs doublet model can enhance the decay width 
cind resolve the dimuon charge asymmetry puzzle, consistent with all the experimental 
constraints. The enhancement of requires large Higgs couplings to the first and second 
generations of quarks, which may help explain the recent 3.2a Wjj excess observed at CDF. 

This work is partially supported by DE-FG02-95ER40917 (B.D.), STDF grant 437 and 
also ICTP Proj. 30 (S.K.), Excellent Research Projects of National Taiwan University under 
grant number NTU-98R0526 (YM) and DE-FG02-91ER40626 (Q.S.). 
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